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Abstract 

We identify the strong coupling fishnet diagram with a certain Ising spin configuration in the lightcone 
worldsheet description of planar Tr$ 3 field theory. Then, using a mean field formalism, we take the 
remaining planar diagrams into account in an average way. Since the fishnet spin configuration is regular 
but non-uniform, we introduce two mean fields (j>, <t> where the fishnet diagram is the case <fi = 1, 4>' = 0. 
For general values of these fields, the system is then approximated as a light-cone quantized string with 
a field dependent effective string tension T e fr (</>, <f>'). We also calculate the worldsheet energy density 
£(4>,4>'), and find the field values that minimize it in the presence of a transverse space infra-red cutoff 
e > 0. The criterion for string formation is that the tension in this minimum energy state remains non- 
zero as e — > 0. In the most simple-minded implementation of the mean field method, which neglects all 
short range correlations of the Ising spins, we find, in this limit, that the tension vanishes for weak and 
moderate coupling, but for very large coupling does indeed stay non-zero. However, a more elaborate 
treatment, taking temporal correlations into account (but still neglecting spatial correlations), removes 
this "phase transition" and the string tension of the minimum energy state vanishes for all values of the 
coupling when e —* 0. Our mean field analysis thus suggests that the "fishnet phase" of Tr$ 3 theory is 
unstable, and there is no string formation for any value of the coupling. This is probably a reasonable 
outcome given the instability of the underlying theory. It is encouraging for our method, that an approach 
designed for a string description can predict, where appropriate, the absence of string formation within 
an intuitive and simple approximation. 
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1 Introduction 



Ever since Maldacena proposed that IIB superstring theory on an AdSsxS5 background is equivalent to 
supersymmetric Yang-Mills field theory with extended M = 4 supersymmetry [1], the status of string theory 
in physics has dramatically changed, from a speculative "theory of everything" with no experimental sup- 
port, to a potentially powerful tool for analyzing the nonperturbative physics of experimentally established 
quantum field theories. The most exciting application of this new string theory technique is, without doubt, 
the confinement problem of QCD. 

While the vast bulk of the literature on this exciting new program of research is dedicated to a "top- 
down" approach, seeking to recover field theoretic physics from the string formulation [2,3], the only work 
toward a "bottom-up" approach, seeking a direct construction of a string worldsheet formalism that sums 
the planar diagrams of quantum field theory, has been the series of articles [4-6] . These three articles are 
foundational for this more practical approach: the first sets up a worldsheet formalism for scalar <E> 3 theory, 
the second for Yang-Mills theory, and the third for the whole range of interesting supersymmetric Yang-Mills 
theories, including TV = 1,2,4 extended supersymmetry. 

The next task is to develop an effective and tractable framework for extracting nonperturbative physics 
from this new worldsheet formalism. The first steps toward this goal were taken in [7, 8] . These two articles 
develop alternative mean field methods on the worldsheet for approximating the sum over planar diagrams. 
The first of these replaces the Ising spins in the formalism by a homogeneous mean field, whereas the second 
makes a mean field replacement of scalar bilinears of the target space worldsheet fields. We continue the 
development of both approaches in the present article by considering a non-homogeneous mean field which 
can be used to describe the fishnet diagram. In this way we can assess the role of this special diagram in 
the context of an average treatment of the sum of all planar diagrams. 

In Section 2 we briefly review the worldsheet formalism for the Tr$ 3 theory. We include a new improved 
way of incorporating the mass of the scalar field. The review closes with Eq. (15) which explicitly displays 
the worldsheet system that sums the planar graphs of this theory. The rest of the article uses mean field 
techniques to approximately solve this system. Here we want to draw attention to the essential features of 
this system, including an explanation of the cutoffs employed to completely specify the dynamics. 

The worldsheet we discuss is based on light-cone parameters, an imaginary time r = ix + = i(t + z)/V2 
in the range < r < T and a worldsheet spatial coordinate < a < P + chosen so that the P + density is 
uniform. Both of these parameters are discretized: r — ka and a = Im with T = Na and P + = Mm. The 
limit of a continuous worldsheet is equivalent to the double limit M,N — ► oo with N/M = (T/ P + )(m/a) 
fixed. In this limiting process the quantity m/a which has the dimensions of force or tension can have any 
value. Thus we actually have a family of cutoff theories labeled by this ratio. The so-called DLCQ cutoff 
corresponds to a — > at fixed m. The continuum physics must be independent of this ratio. 

Qualitatively the worldsheet formalism maps every planar diagram to a worldsheet with several internal 
boundaries each at fixed values of a. The worldsheet target space field q(a, r) satisfies Dirichlct conditions 
q = ki on the ith boundary, and each ki is integrated (these are the loop momenta of the multi-loop 
diagram). On the discretized worldsheet the boundaries are also discretized, occupying a certain number 
of temporal links. When a boundary occupies a temporal link — 1), Dirichlct conditions require 

the insertion of a factor 8{q{ — ql~ ) in the path integral. To keep track of the occurrence of boundaries, a 
two valued Ising-like variable P- — 0, 1 is assigned to each site. The value 1 means the site is crossed by a 
boundary and means it is not. Thus each consecutive pair of sites crossed by a boundary comes with the 
factor 

in the path integral. If cither of the P's here is zero the factor reduces to unity. Certainly for numerical 
work and also for conceptual clarity, we choose to use a Gaussian approximation to the delta function: 
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If wc keep e finite, we see by comparing the two right sides at q\ — q\~ x = 0, that (2ira/me) d / 2 = L d 
where L is the size of the transverse space. That is we have an infra-red cutoff L in transverse space 
L 2 = 2ira/me. Exploring the connection between field and string theory in the presence of an infra-red 
cutoff is not unprecedented. Indeed, the work by Berenstein, Maldacena and Nastase [9] which explores 
the field/string connection in the original SUSY Yang-Mills/AdS context defines the gauge theory on S3 
rather than i? 3 . Similarly, in the present work we study the sum of planar diagrams with e finite. We must 
remember though that the "field theory" we are thereby solving is some compactified version of the original 
$ 3 theory. The worldsheet formalism for a more conventional compactification (on a transverse torus) of 
Tr<f> 3 theory is sketched in Appendix D. 
Using (2) leads to the factor 

-j expj— ^J-W-^ (3) 

for each pair of sites crossed by a boundary. One can immediately anticipate that the e — > limit will be very 
delicate in approximate calculations. For example the mean field approximation we shall employ replaces 
P- by a slowly varying field <f> taking values between and 1. The q integral then behaves as e d / 2 whereas 
the prefactor is softened and the net factor becomes e^ 1- * 2 ^/ 2 , which strongly suppresses the contribution 
as e — * except when </> = 1. We can mitigate this delicacy somewhat by arranging the e factors as in 
(15). This makes no difference in the exact formula but makes a big difference in the e — > limit of a 
mean field approximation. Even if the P's are less than 1 the singular e dependence cancels between the q 
and b, c integrals. This helps define a less singular mean field approximation, but the need for such a finely 
tuned set up of the approximation nonetheless casts doubt on its reliability at small e. We can be more 
confident about our conclusions for e relatively large, say of 0(1). It is nonetheless important to confirm (or 
refute) our conclusions with more exact calculations such as Monte Carlo and block spin renormalization 
group studies. In any case we think our mean field studies give an important qualitative description of the 
worldsheet dynamics that will be crucial in interpreting these more accurate numerical methods. 

In Section 3 we give a brief description of the fishnet diagram and its representation in terms of a 
specific Ising spin pattern on the BT worldsheet. This pattern motivates the choice of mean fields we use to 
illuminate the role of the fishnet in the complete sum of planar graphs. Sections 4 and 5 explain the results 
of our two approaches to the mean field approximation based on mean fields (j? k — (P^) and ipj. = (P^Pfe -1 ) 
respectively. Section 6 gives our concluding remarks. Appendices A, B and C contain the technical details 
of our calculations. 



2 Worldsheet for Summing <£ 3 Planar Diagrams 

We start our review of [4] by recalling the worldsheet representation of the free gluon propagator (we shall 
refer to the 4> field quanta as "gluons" even though their spin is 0). Consider a gluon carrying M units of 
p + = Mm and and propagating N time steps. The transverse momentum p = q^ — q , where the worldsheet 
target space fields q 3 k , k = 0, . . . , M and j labeling the time slice, live on the sites of the worldsheet lattice. 
q\ satisfies Dirichlet boundary conditions <7q = <7 , q ] M — q M . We also must introduce Grassmann ghost 
variables &i,c{ on each lattice site to provide crucial factors of M when we represent a single gluon as M 
bits. D = d + 2 will denote the space-time dimension, which we keep general for our formal presentation. 
But in our numerical work we take the case d = 2, four dimensional space-time. There are d components of 
each q and d/2 pairs b, c of ghost fields. With this notation the action is 

S = S g + S q (4) 

M-l 

j »=o 
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s B = --E 



M-2 



(6) 



i=l 

Then the master formula for the worldsheet formalism is [4] 

ex P {-g!} = expj-^^o)!} = /n|^ ee / DcDbDq (7) 

In (7) we have imposed Dirichlct boundary conditions. However, to dynamically implement momentum 
conservation in the path integral we integrate over q 3 M independently, retain the Dirichlet boundary condition 
at i = 0, where we impose q 3 = q , but insert momentum conserving delta functions at the other end i = M. 
It is convenient but not necessary to set q Q = 0. 

We define the dimensionless coupling 

Note that in D space-time dimensions we have all together d/2 sets of b, c ghosts, denoted by bold-faced 
letters when referring to all such sets. But one of these sets, identified with non-bold type, is singled out to 
handle the special factors occurring at interaction points. When needed, we denote the remaining (D — 4)/2 
sets by checking their symbols b, c. 
We associate factors 

V 3 0i ee flexp{-^(6?"_ 1 ^_ 1 + 6j +lC ? +1 )} (9) 

K = ge X p{-^(i4ti-f4 +1 )(4t 1 i-4 +1 )} (10) 

with the fission and fusion vertices respectively, the different forms due to the asymmetric way the ghost 
insertions assign the required l/p + factors to the vertices. 

Next we recall the formula that systematically sums over all the planar diagrams on the lattice. The 
worldsheet for the general planar diagram has an arbitrary number of vertical solid lines marking the location 
of the internal boundaries corresponding to loops. Each interior link j, j — 1 of a solid line at spatial location 
k requires a factor of — q^ 1 )- To supply such factors, assign an Ising spin s 3 k = ±1 to each site of the 
lattice. We assign +1 if the site (fc, j) is crossed by a vertical solid line, —1 otherwise. We also use the spin 
up projector P 3 = (1 + s 3 k )/2. 

At the endpoints of each solid line we have to supply the vertex insertions Vo, Vo- These factors occur 
when there is a spin flip, s 3 k = —s 3 ^ 1 . In the foundational papers we exploited the properties of Grassmann 
integrals to take these factors into the exponent where they were multiplied by bilinears in the ghosts. This 
was particularly convenient for the vertices of gauge theories, which had momentum and spin dependent 
prefactors. However, for the scalar theory treated in this article, this is not necessary, because the prefactor 
is simply g and the ghost insertions are already in exponential form. In the general formula we need only 
multiply the be term in the exponent of the vertex insertion by P 3 P 3+ (1 — P 3 ~ ) or P 3 P 3 ~ (1 — P 3+ ) 
respectively. Note that each internal loop occupies at least two time steps. 

In this article we implement the Dirichlct conditions on boundaries using the Gaussian representation 
of the delta function (see the last line of (2)). We keep e finite until the end of the calculation. Using this 
device, our formula for the sum of planar diagrams becomes: 

Tfi = lim E fDcDbDqcxp\\ngY,( P ' P ' +1 ^- p r 1 ) + ^Pr 1 ^-Pr))\ 

sl=±l [ i,3 J 



exp 



3 




y i,j 

ex p I £ YM+i 4m + , w - w - ^+i) + £ £a - + /v J> (n) 



The first exponent in this formula supplies a factor of g whenever a boundary is created or destroyed. The 
second exponent includes the action S q for the free propagator together with the exponent in the Gaussian 
representation of the delta function that enforces Dirichlet boundary conditions on the solid lines. The 
contents of this exponent go to the discretized action for the light-cone quantized string, if the quantity 
a 2 P- P{~ X /m 2 e is replaced by 1/Tq, with T the string rest tension. The third exponent incorporates the 
e dependent prefactor for that representation of the delta function as a term in the ghost Lagrangian. 
The remaining exponents contain S g together with strategically placed spin projectors that arrange the 
proper boundary conditions on the Grassmann variables and supply appropriate 1 /p + factors needed at the 
beginning or end of solid lines. 

We remark that, when Dirichlet conditions are imposed at i = 0, M, the expression (11) sums all the 
planar multi-loop corrections to the gluon propagator. The evolving system is therefore in the adjoint 
representation of the color group. That is, we have tacitly assumed that the only solid lines initially and 
finally are those at the boundaries of the strip. More general initial and final states are described by allowing 
more solid lines initially and finally. When the system is in a color singlet state (a "glueball"), we must 
include diagrams in which the outer boundaries are identified, i.e. the diagrams should really be drawn on 
a cylinder, not a strip. In this case, q(p + ) = q(0) +p, the variable q(<r) is strictly periodic only in the state 
of zero total transverse momentum. 

When we regard the worldsheet as a cylinder, we realize that the Ising spins can be in configurations 
that have no correspondence to the original Feynman diagrams. These involve time slices in which all spins 
are down, signifying the complete absence of solid lines through the slice. Several consecutive such slices 
represent the evolution of a no gluon or "nothing" state. This is the closest thing the formalism gets to a 
closed string state in the non-interacting theory. The worldsheet formalism associates this state with the 
ambiguous 0/0: both the be and q systems have a zero mode on a nothing time slice, the ghost zero mode 
is responsible for the in the numerator and the q zero mode for the zero in the denominator. In the exact 
formula this nothing state decouples from an initial state with gluons because the ghost insertion associated 
with a terminating boundary supplies a second zero factor making the contribution of the transition time 
slice 2 /0 = 0. So the nothing state will be removed from all amplitudes that have at least one gluon 
in the initial state. However if one is considering the generic state that evolves without regard to initial 
conditions, the nothing state will be included. In particular, the mean field approximation employed in this 
article replaces the Ising spins by average spins, the approximate amplitudes will not involve definite spin 
configurations, and for non-zero mean fields the zero modes are removed. The nothing state will be mixed 
in with multi-gluon states. For moderate to strong coupling where the typical time slice has many up spins 
this admixture should be innocuous. However, for weak coupling when up spins are very rare, the precise 
connection with straight perturbation theory will be obscured because the mean field method will treat the 
precise decoupling of the nothing state in a rough average way. So applying the mean field approximation at 
nonzero coupling and then taking the coupling to zero will most likely lead to the description of the nothing- 
state rather than a state with a small number of gluons. 

After that digression on "nothing" we continue with our discussion of the world sheet path integral. We 
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can simplify the first exponent with the replacement 



exp 



{in^(^^ +1 + 1 - ii': I'-: 'r ] - H (5) (1 - s * srl)/2 . (12) 



This is valid because the left side agrees with the right side except for spin patterns that include one or 
more isolated up spin (one preceded and followed by a down spin) . But inspection shows that the integrand 
is independent of the b 1 c on the site of each such up spin, so the ghost integral gives zero. Thus only spin 
patterns in which each up spin is preceded and/or followed by an up spin contribute to the worldshcet 
integral. This also explains why the original formula vanishes at zero coupling if there is one or more spin 
flip in the spin pattern. 

The formula (11) assumes the $ field is massless. We now present a new worldsheet local procedure for 
incorporating a mass fi for $ that is superior to the one offered in [4]. Consider the identity 



1 dcidbi \ a 
— ?. — exp < — 
2ir I m 



n 



blCi : &M-1CM-1 
1 



M-2 



+ —7— + ^ (6 i+1 - bi)(c 




(13) 



We have already exploited two special cases of this formula: a = (3 = 1 gives one time slice of the worldsheet 
ghost path integral for the massless propagator, and the case a = 00, (3 = 1 provides the mechanism for 
supplying l/p + factors in the worldsheet path integral for interacting Feynman diagrams. 
A finite mass [i in the free propagator corresponds to a factor 



2p+ J V dMm) 

where the right side is a valid discretization T = Na,p + = Mm of the left side, going to the latter in 
the limit M, TV — > 00. Thus we must provide a factor (1 — p 2 a/dmM) d ^ 2 on each time slice. This will be 
achieved with d/2 copies of (13) if we choose a + (i — 2 = —p 2 a/dm and multiply the left side of (13) by 
(a(3) d / 2 = cxp{(d/2) In a/3}. For example, with /3 = 1, the choice a = 1 — p 2 a/dm does the job. Define 
the parameter p = (1/a) — 1 = p 2 a/(dm — p 2 a). Then the worldsheet action for the massive propagator 
is obtained by adding the term — p&jcj to the free ghost action (6) and also multiplying the path integral 
by a Nd / 2 = (1 + p)- Nd ' 2 . This latter factor can be associated with the boundary of the world sheet, since 
its exponent is proportional to the length of the boundary. It contributes like a "boundary cosmological 
constant" . We should include one such factor for every internal boundary and the square root of the factor 
for each external boundary. (When we draw our graphs on cylinders every boundary will be internal.) Since 
boundaries occur whenever Pi = 1, for a general diagram we can supply this factor by including the term 
(d/2) ^1 m tne exponent of the general worldsheet integrand. 

In the worldsheet path integral for a general diagram, the added term pb^c^ should appear on a site with 
Pi = when it is immediately to the right of a solid line, that is when Pj._ l = 1 as well. Furthermore, it 
cannot occur on such a site if it immediately follows the beginning of solid line, since the coefficient of be 
on such a site must be zero (to properly include l/p + factors). Thus we must have P^Zi = 1- The upshot 

is that the term should be multiplied by the factor (1 — Pt^Pk-i^l-i- Because of these restrictions, we 
might wonder whether the contribution of the boundary factor should share them. However, we think not. 
First of all, two adjacent solid lines bound a gluon carrying a single p + unit, i.e. M = 1. According to (14) 
an M = 1 gluon should have a factor a d l 2 for each time slice. This is precisely supplied by the boundary 
factor. (Since an M = 1 gluon contains no spin down sites, there is no pbc term.) So we should include the 
boundary factor for every solid line, including those adjacent to each other. Next, what about the absence 
of the pbc term at the beginning of each solid line; should we also omit the factor a d / 2 for that one time 
slice? In other words, should a boundary of N time steps yield the factor a d ^ N ~ v> / 2 or a dN / 2 7 Actually 
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this is a matter of taste because the discrepancy between the two choices can be absorbed in a redefinition 
of the coupling constant. In order to preserve the interpretation of these factors as boundary cosmological 
constant effects we choose the second alternative a dN ^ 2 . 

After all of these considerations we propose that the worldsheet path integral for the sum of planar 
diagrams in the massive theory be given by 

j j—i 



Bij 



Tfi = lim £ J DcDbDqexp J mg^ 1 ^ - d - In (1 + p) £ P{ 



exp<J-^ «,,fr!r!+C,(>! +J - t^(c^, - c/ ) - D, : , l>> ),r ; ' +l - <■> ) 

Mi = ^/r/y : • i'; r-! r; > • Q f?')(i? +] i-itiH pO pM ,/V : (16) 



e 



a - pi) (pLitiSi 1 - p i+i)+ pun-i^ -pi~i)+ pt'pt'pL) a?) 
ca = a-pi)(i-n +l ) (is) 
Da = a-pi)(i-pi + jpr'pt 2 (19) 

In the rest of the article we apply mean field methods to this version of the worldsheet system. It is helpful 
to note that once /x 2 ^ 0, one can scale p to p = 1 by fixing an appropriate value of m/a. Then the massless 
limit would be regained by taking g — > oo. 



3 The Fishnet Spin Pattern 

We would like to extend the mean field formalism of [7, 8] to allow for regular spin patterns that are not 
necessarily completely uniform. In spin models, such a generalization is important for describing the anti- 
ferromagnetic phase. Similarly in the sum of planar diagrams, the fishnet diagrams first identified by Nielsen 
and Olesen and of Sakita and Virasoro [10] are an important subclass of diagrams that naturally contain 
a stringy interpretation. Based on these ideas the program of using a "strong coupling fishnet" diagram 
as the scaffolding for the construction of a worldsheet representation of the sum of the planar diagrams of 
quantum field theory was proposed in 1977 [11]. This work made use of the mixed (x + ,p + ,p) representation 
[12, 13] of each propagator, and further employed a discretization of both p + = Im and ix + = ka with k, I 
positive integers [14]. In the context of this discretization, the formal strong coupling limit singles out those 
diagrams in which every propagator spans a single time unit a and carries one or two units m of p + . Among 
such diagrams, the fishnet diagram is further restricted by the requirements (1) that it be planar, and (2) 
that no pair of propagators connect the same pair of vertices. These properties together single out a more 
or less unique diagram, the strong coupling fishnet. For Tr$ 3 field theory, this diagram, calculated in [15] 
is shown on the left in Fig. 1. It was shown in [11,15] that the diagram describes the propagator of the 
lightcone quantized string [13], discretized as in [14]. 

Of course there are many other diagrams that are important at strong coupling: the fishnet diagram was 
to be regarded as a worldsheet template for the inclusion, in principle, of all the rest of the planar diagrams. 
The proposal amounted to a nontrivial reorganization of the summation of planar diagrams in a fashion 
in which a string worldsheet description occupied center stage. However, in the references cited above, no 
approach to the inclusion of "the rest of the planar diagrams" , beyond a brute force calculation, was offered. 

The worldsheet formalism [4] , reviewed in the previous section, offers powerful techniques for an approx- 
imate evaluation of the contribution of all planar diagrams. We begin by identifying which pattern of Ising 
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Figure 1: The spin pattern on the right corresponds to the Tr$ 3 fishnet diagram on the left. Both diagrams 
are drawn on a cylinder, with opposite sides identified. 



spins corresponds to the fishnet diagram just described. We have drawn the corresponding Bardakci-Thorn 
(BT) worldsheet on the right in Fig. 1. The Ising spin at a given site is +1 or —1 depending on whether the 
site is crossed by a solid or dotted line respectively. Following the sites for a fixed spatial value, we find the 

repeated pattern — h + H h + H !- + + •••. The corresponding pattern at the neighboring spatial sites 

is + H h + H h + H h • • •• This has an anti-ferromagnetic character, though not quite the strictly 

alternating pattern of the classic Ising anti-ferromagnet. 

Our calculations are based on the mean field method, developed for this system in [7, 8]. The first of these 
articles assigned the mean field <f) 3 k — {Pi), and studied the energy for a completely homogeneous pattern 
4>° k = 4>. In Section 4 we generalize this work to the fishnet spin pattern by introducing two mean fields <f), <f)' , 
with <j) the average spin at the sites marked by + in the previous paragraph and with <jjf the average spin 
at the sites marked by — . Note that with <f> ^ <f>' the translational symmetry of the lattice is broken by the 
assignment. But the homogeneous spin pattern can still be accessed as the special case 4>' — <j>. We shall 
calculate the bulk energy per site £((/>, 4>') of the system, as a function of <f> and 4>' . The fishnet diagram 
corresponds precisely to the case <f) = 1 and <f>' = 0. Configurations away from these values take into account 
all other diagrams in an average way. Our work therefore provides a first look at how the complete sum of 
planar diagrams alters the string interpretation given by the fishnet. 

The second mean field article [8] applied the mean field method to the matter fields rather than to the 
Ising spins. However, this approach also leads to a mean field description of the spins in which the mean 
field ipi = (Pi 1 P'l), and is associated with temporal links rather than sites. For a link joining two + 
spins, the field has the value 1, whereas it is zero for all other possibilities. For the fishnet pattern we then 
find the link field values 01100110011 • • • alternating with 10011001100 • • • on spatial slices. Although this 
second approach is in some ways superior to the first, its application to the full-fledged theory is considerably 
more complicated than the first and won't be attempted here. Instead, we explore the role of the spin order 
parameter tf? k = (PlP^ 1 ) singled out by this second approach. In section 5 we simply repeat the calculations 
of Section 4 using this alternate choice for mean field. 
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4 Mean Field Treatment of the Ising Spins 



In this section we work directly with the Ising spin system and approximate the P/'s in certain terms of 
the action by a mean field as was done in [7]. We do this by first setting up the standard effective action 
formalism, by adding a source term ^2 k - n k P k for the P variables to the action. Then we write the path 

integral in the presence of k as e w ^ and define the effective field <j? k = (Pl) K — dW/dn k . The effective 
action is then obtained by Legendre transformation 

A{4>) ee W ( K )-J244 (20) 

kj 



dA 



(21) 



Then it follows that, when the source is removed by setting it to zero, the spin expectation <j> is a stationary 
point of the effective action. 

So far everything is exact. The mean field approximation comes into play in attempting an approximate 
calculation of W(k). A straight perturbation theory would treat the terms coupling the matter fields and 
the spin fields as a perturbation. The way we implement the mean field approximation is to improve on 
perturbation theory by replacing the spin dependent coefficients of the matter terms in the action by their 
expectation values, regarded as fixed for the purposes of the spin sum. This decouples the spin sum from 
the matter integrals and makes it tractable. This decoupled spin sum then yields an approximation to the 
k dependence of W, which in turn determines the expectations (P), (PP 1 ), (PP'P"), which are used as the 
zeroth approximations to the coefficients of the matter terms. One then looks for stationary points of the 
approximate total effective action. This zeroth order mean field approximation can then be improved by 
treating the differences between the actual spin dependent coefficients and their assumed expectation values 
as a perturbation. 

The extent of the approximation depends on how much of the P dependence is replaced by its expectation 
value. The simplest option is to make this replacement for all terms in the action (excluding the source term 
itself of course). But we shall also (as in [7]) consider replacing only the spatially coupled P's and those 
multiplying matter fields by their expectations, since the remaining terms can be explicitly included in the 
spin sums 4 . The first option neglects all short range correlations between spins, whereas the second takes 
at least some of the temporal correlations into account. 

In all cases the matter (q and be) integrals are evaluated in the presence of the mean fields <j). These 
integrals are done in Appendix C for the fishnet pattern involving the two mean fields <p, <fi' . In our version of 
the mean field approximation the monomials of P's in the matter action are replaced by their expectations. 
The fishnet spin pattern requires only five such quantities: (++), ( — h), (+++), (-++), (h — +), all defined in 
Appendix A. The result of the q integration (Eq. 177) is repeated here for convenience. 



(2Trm/a)- dMN/2 e w «= JJ 

KM/2 



{a 2 - /3 2)W-2 sinh ^+ sinhA^ 



l i- d / 2 



sinh2^ sinh2£i 



In view of the N = T/a dependence in this formula it is clear that the £± represent excitation energies of 
the system: 

AP± = -£±. (22) 

4 In the standard mean field treatment of the Ising model every spin in the action is replaced with a mean field. But replacing 
say only the spatially coupled spins by mean fields leads to an improved approximation to the critical temperature. In the two 
dimensional Ising model the usual mean field treatment gives T™^ f» 1.76T^ xact . In contrast replacing only the spatial spins 
by mean fields leads to the improved estimate r pP artlal ~ 1.55T£ xact . As is well known, neither version does a good job with 
critical exponents. 
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Finite energy excitations correspond to = 0(1/M), which entails y/Ri sinhK; = 0(1/M) or y/Ri coshK; 
1 + 0(1/M 2 ). Since ui[ ~ (ln/M) for I < M we read off in this limit 



^R~i sinh ^ V7 



2 + a 



(23) 



M V 2a + a 2 - (3 2 ' 

Since Amass? = 2MmAE = 2nT cS l we read off the effective string tension in terms of the background fields: 



m 

Jeff = — 



2 + a 



m 



(++> + (-+>+4e 



a V 2a + a 2 - (3 2 a y 2«++) + (-+) + (++)(-+) /e) ' 
where a = ((++) + (-+})/2e and /3 = ((++) - (-+))/2e. 



(24) 




Figure 2: The contribution of g plus ghost fields to the bulk energy. In the left graph we have taken e — > 
holding <j)' = ef fixed. The right graph is plotted for fixed e = 0.001. Note that the effect of finite e is 
concentrated near 4> = 0. 

For fixed <p, <f>', both a and (3 are of 0(1/ e) as e — > 0. With no relation between them, it is easily seen that 
the tension is of 0(y/e) and hence vanishes. If, however the limit can be taken in a way that (a 2 — /3 2 )/a is 
finite, the tension will also be finite. This can be achieved either with both (++), ( — h) = 0(e) or (++) = O(l) 
while (-+) = O(e). In the second case, writing = ef and ( — h) = efG(<j>), we find 



a^2(l + fG{4>)) 



(25) 



For generic fields the mean field approximation amounts to a string description of the system, with a field 
dependent tension. But only if fG(<fi) is finite for the field values that minimize the energy, will the system 
behave like an actual physical string. 

The matter contribution to the ground state bulk energy is (see Eq. 184): 



■ ln(l + p) 



£ m (*,*0 = ~(€i+6)+2]jf E Ki+d+ln(a 2 -/3 2 )] + ^ 

KM/2 

Since the explicit expression is quite complicated, we display this function graphically. In Fig. 2 we show the 
matter energy as a function of <f> for a sample of <j>' values. In these graphs we are assuming the factorization 
(++) = (j) 2 , (-+) = 4>4> ', (+++) = (f> 3 , ( — h+) = (h — h) = 4> 2 4>'- The effect of the infra-red cutoff e is seen by 
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Figure 3: The contribution of the spin sum to the bulk energy per site using (30), which neglects all short 
range correlations. 



comparing the case with e = 0, / = <fi' /e fixed, to the case with e = 0.001, <j>' fixed. It clearly has the most 
dramatic effect in the region of small <f>. 

As already mentioned, we present two treatments of the Ising spin contribution to the energy. In the 
first (simplest but most drastic) option all spin dependence in the action is replaced by its corresponding 
expectation value, and the spin sum reduces to a trivial independent sum on each site: 

Y[ J2 e< p /=exp|^ln(l + e<)l. (26) 

kj P J k =0,l {. k i ) 

This version completely neglects short range correlations. Since this represents the entire k dependence of 
the approximate W the relation of <fi to k is independent on each site: 



1 + e< • 7 



Because the spin sum is independent on each site all expectations of the monomials of P's appearing in the 
action factorize into the corresponding monomials of (j>. 
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We denote by A s the spin part of the effective action, including the — K<f> term of the Legendre transform 
and also the terms involving the coupling constant: 



■As = ^[-44+ln(l + e<)]+21n^((P^(l-^ +1 )) (28) 

kj kj 

]T [-4 ln# - (1 - 4) ln(l - 4)] + In? 2 ]T 4(1 - 0£ +1 ) (29) 



-MAT j 2 [0 In + (1 - 0) ln(l - <j>)} + ^ [0' In 0' + (1 - 0') ln(l - 0')] 
-1 (30 + 0' - 200' - 20 2 ) In? 2 J = -MN£ s (0, 0'), 



(30) 



where the last line specializes to the fishnet spin pattern. Since the pattern is regular in the temporal 
direction the coefficient of — N is just the energy in units of 1/a and the coefficient of —MN is the energy 
per site £ s . 

We display the spin energy per site for various couplings in Fig. 3. These graphs display three distinct 
phases for the spin system. At weak coupling the 0' = curve indicates a degenerate ground state with 
~ 0, 1. These values correspond to the Ising spins aligned with values —1, +1 respectively, a ferromagnetic 
phase. The curves with 0' constrained to 0,0.1,0.2 approximate the 0' = curve more closely as 0' — > 0. 
At intermediate coupling the minimum of the 0' = curve is fixed at = 1/2, corresponding to (sj.) = 
for the Ising spins. This is a disordered phase. The constrained 0' curves give more information about 
this phase. At g 2 = 1 they all show a minimum for = 1/2. For g 2 > 1 the minima of these curves are 
above 1/2, approaching 1/2 as 0' increases. In other words the sites are anti-correlated with the 0' sites. 
This anti-correlation will tend to reduce the paramagnetic susceptibility. For g 2 < 1 (not shown) the curves 
would have minima less than 1/2 that all lie on the 0' = curves. In other words the sites are positively 
correlated with the 0' sites, and this will tend to increase the paramagnetic susceptibility. In summary, 
although the phase is not ordered, it shows a paramagnetic susceptibility that increases as g 2 decreases. 
Finally for strong coupling the 0' = curve no longer contains the minimum energy. The minimum energy 
is on the curve with 0' = at ~ 0.7, an anti-ferromagnetic phase. This phase structure is not surprising 
because the spin energy is essentially the mean field free energy of an Ising model, and the phases we have 
described are well understood properties of the Ising model provided its dimensionality d > 1. The mean 
field description fails to predict the critical exponents for d = 2. 

However, the spin sums we have actually done correspond to Ising couplings in one dimension only (the 
temporal direction), and it is well-known that the Id Ising model displays no phase transitions: there is 
only a single disordered phase, albeit with paramagnetic properties. In our worldsheet system there are 
some direct spatial couplings of the spins, and there are also indirect spatial couplings induced through the 
spin-matter coupling terms. When these spatial couplings are brought into play, a phase transition is no 
longer ruled out, but it still behooves us to do a better job on the one-dimensional spin sums that removes 
this artificial phase transition. For that reason we now turn to the more sophisticated treatment of the 
spin sums sketched in Appendix A. In this better approach the spins multiplying In g are not replaced by 
mean fields but rather incorporated exactly in the spin sum. However this sum is tractable only after one 
specializes to a regular spin pattern, here the fishnet spin pattern. Then the problem is solved by finding the 
eigenvalues of the 2x2 transfer matrix and choosing the largest eigenvalue. In the notation of Appendix A 
the spin energy per site now becomes 

£ s = ~lnt%(K, k') + ^k+^K^+^k'. (31) 

Again we display the spin energy in this case graphically in Fig. 4. We see that the indications of a phase 
transition seen in Fig. 3 have now disappeared, since we now have an essentially exact treatment of the Id 
Ising model reflected in the spin sum. Now the minimum energy is always on the 0' — curve and located 
at = 1/2, the disordered phase. But the 0' ^ curves still show that this phase is paramagnetic: for 



11 



Figure 4: The contribution of the spin sum to the bulk energy, using (31) which includes coupling de- 
pendence, and hence some temporal correlations, in the spin sum. Note the disappearance of the phase 
transitions. 

g 2 < 1 the susceptibility is greater than in the disordered phase and it is less for g 2 > 1. We shall see that 
these qualitative features are maintained when we add in the energy of the matter sector. The <f>' = curve 
follows the fate of the fishnet diagram, corresponding to the endpoint <j> = 1, when the rest of the planar 
diagrams are included. In the following analysis we shall always display plots for the total energy per site in 
this second more accurate treatment of the spin sum. 

We now come to the total energy. The details of the matter contribution to the energy are contained in 
Appendix C. For ease of comparison we shall display graphs for the same parameters we chose in the spin 
energy. We start by displaying the total energy curves for e = 1 so that the infra-red cutoff in transverse 
space is L ~ 0{y/a/m). These results are shown in Fig. 5. We first note that as with the spin energy, it 
is still true that the lowest energy solution always falls on the <f> = <fi' curve, and there is no indication of a 
phase transition. However, the location of the minimum now depends dramatically on the value of g 2 . It 
evidently approaches 4> = for weak coupling and (j> = 1/2 for strong coupling. Looking at the <f>' = curves 
we see first that the fishnet diagram <p = 1 on this curve predicts positive energy for weak and intermediate 
coupling. By g 2 = 100 this diagram predicts negative energy. But in all cases the minimum of the </>' = 
curve is negative but higher in energy than the uniform 0' = <f> case. 

Next we display graphs with e = 0.00001 and the same range of couplings. We see from Fig 6 that 
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Figure 5: The total energy for e = 1. The location of the minimum depends dramatically on the coupling. 

qualitatively the energy curves resemble more closely those of the spin energy alone. The minimum energy 
again lies on the = curve and is located at (j) = 1/2 for all couplings. And again we see the same 
paramagnetic tendencies in the fixed <f)' curves for strong and weak coupling. The string tension for this 
lowest energy state can be found explicitly. Since (j) = <t>' = 1/2 at this point, it follows that n = so 
(++) = (-+)'= <f>/t+(0) = 1/(2 + 2.9). Thus 

T oS = ^y/2e(l+g), (32) 

which vanishes for e — > at fixed g. Vanishing tension implies among other things a gapless energy spectrum. 

In Fig. 7 we display the energy curves for a massive scalar field. It is interesting that, for weak coupling, 
there are indications of a ferromagnetic behavior. (We were able to do calculations at zero mass because of 
our cutoffs.) Notice though that for g 2 > 1 the massive case is very similar to the massless case. In particular, 
at strong coupling we see a continuous approach to the disordered phase 0=1/2 and no anti-ferromagnetic 
behavior. Notice that although there are quantitative differences, the energy curves for the massive case at 
e = .00001 are qualitatively similar to those of the massless case at e = 1. This is perhaps reasonable because 
both p = 1 and e = 1 set an infra-red cutoff of the same order of magnitude L 2 = 0(a/m). 

Let us now consider what these calculations tell us about the role of the fishnet diagram. We can follow 
the fate of the fishnet diagram in two stages. Starting at (j)' = 0, <fi — 1 we first let <fr relax at fixed (j)' = to the 
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Figure 6: The total energy for e = 0.00001, corresponding to a very small infra-red cutoff (very large 
transverse space). The location of the absolute minimum is at = 0.5 for all coupling, but the paramagnetic 
properties of the spins show a clear dependence on the coupling. 

minimum energy at cf>' = 0. Since ( — h) = 0, the string tension along this curve is {m/a)\J (1 + 4e/ (++))/2, 
which is simply m/(a^/2) as e — > 0. Then one can begin to relax cj)' putting 4>' = ef. Then one can smoothly 
send e to zero at fixed /, and study how the physics evolves as / increases. The field dependence of the 
energy is shown in Fig. 8. One should keep in mind that because of our limiting procedure all these curves 
are infinitesimally near the 4>' — curves of the previous graphs. In Fig. 9 we show the string tension 
corresponding to the minimal point on each fixed / curve as a function of /. Next to it we display the 
corresponding minimal energy as a function of /. We see that both string tension and energy generally 
decrease as / increases. Note, however that a small energy barrier between / = and / > 10 is visible on 
the g 1 — .01 curve. Although not visible there are similar tinier barriers at higher coupling. It is clear from 
the previous graphs that the general decreasing trend continues till infinite / where the tension goes to zero 
and the energy reaches its minimum. Our conclusion is that, for Tr<E> 3 theory, the fishnet phase is not stable 
and though there is a stringy description for generic values of the fields, the system relaxes to a ground state 
with no mass gap and zero string tension. In the table we have compiled a list of the relevant parameters 
characterizing the physics of this global minimum for various g, p, e. 

What would we expect to happen in a theory like QCD which ought to confine quarks by forming stringy 
flux tubes? The failure of the Tr$ 3 theory to form string isn't just signified by the monotonic decrease of 
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Figure 7: The total energy for the massive case, p = 1, e = 0.00001. There is evidence for ferromagnetic 
behavior at weak coupling but none for anti-ferromagnetic behavior at strong coupling. Note the qualitative 
similarity of these curves to those of Fig. 5. 

effective tension with /, but rather with the correlation of this behavior with the monotonic decrease of the 
associated system energy. If the system energy had a minimum at finite /, its ground state would support 
a non-vanishing string tension at the corresponding /. This is what we should expect to happen for QCD 
and other confining theories. Moreover, for QCD this should occur for couplings of 0(1). 

5 Mean Field Treatment of the Matter Fields (Large d = D - 2) 

In the previous section we built the mean field approximation on the mean field (f> = (P) representing the 
average value of the Ising spin, following [7]. But another, perhaps superior, approach is to instead focus 
on a mean field defined as the average value of an SO(d) scalar bilinear of the matter fields, as suggested in 
[8]. (Such a description should be accurate in the parametric limit d — > oo. The existence of this parametric 
limit in which the approximation is good is a significant virtue of this approach, although the underlying 
field theories we are describing are only meaningful in sufficiently low dimensions.) In this approach we want 
our mean field to be an order parameter which measures the presence of an effective string tension, so it is 
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.00892 
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-2.3792 
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.0455 


.0041 


.01482 
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5 


-2.0353 


0.4895 


-.4189 


.0360 


.0026 


.0167 
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1 




-4.501 


.4892 


-3.031 


.001 


2.25-10~ 6 


31.62 


10 4 





10" 


5 


-4.574 


.505 


1.8189 


.0124 


3.03-10" 4 


.0284 



Table 1: The worldsheet parameters for Tr$ 3 theory, fixed by minimizing the energy in the mean field 
approximation. 



natural to define 
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me \ 2 

Accordingly, we add a source term to the worldsheet action 



-Y.JI 



me 



kj 



(33) 



(34) 



The path integral in the presence of J is defined to be e w ^ , the field is then xi = dW/dJ 3 k , and the effective 
action is given by 



kj 



(35) 



The mean field approximation is then implemented by establishing the relation between \ an d J as that 
arising from the path integral over matter fields neglecting the terms in the action that couple the matter 
to the Ising spins. In other words we take the J dependence of W to be that of W m defined by 



- / DcDbDqev j £ q{? - ^ £ J{W k q^ f + £ ^ J&cJ 1 

[ k,j k,j k,j J 



(36) 



We then use this approximate path integral to approximate the matter dependence of the terms in the 
action coupling to the Ising spins by their average values. In practice we assume that J (and hence also %) 
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is uniform in k,j. Then by consulting Appendix C, we find that 



W m (J) 



-MN 




m 



J \ - 
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x(J) 



J 



V J 2 + 4e J 



(37) 

(38) 
(39) 
(40) 



Inserting these results into the summand for the spin sum, we find that the spin sum is just a two dimensional 
generalized Ising model with Boltzmann factor e A ^ with 



A(P) = lngJ2 

kj 



d 



+E 

kj 



VJ 2 + 4eJ 



d -A> 

2 kj 



Bl 



k,j 



+ 



kj 



VJ 2 + 4eJ 



2 kj 



kj 



= prn prnpr + a - pim+i + pu) + oo. nmziPLi 
= ^-n)^-H + i)-^=HH + i-pi-n+i 

/P=0.1 e ^ 



and B and D are unchanged. Defining e Ws ^ J ^ — X^p=o l tne total effective action is 

4 w VK m (J) + W S (J) - MNJx(J). 



(41) 

(42) 
(43) 



(44) 



and the value of J (or of x(J)) is determined as a stationary point of A that minimizes the energy per site 
£ = —A/MN. To calculate W s with no further approximation we would need the exact solution of the two 
dimensional spin system defined by (41). This system differs from the soluble two dimensional Ising model 
by linear terms in the spins (a constant magnetic field) and also higher than quadratic terms in the spin. 
For generic J and coupling parameters this model certainly has no analytic solution. A numerical solution 
can be attempted, but we leave this problem for future research. 

Instead we note that, in this setup, the mean field \ ac ts as a source for the spin bilinear (P^P£~ ), and 
so its conjugate J is in this sense a mean field for this bilinear. Thus a natural mean field approach to the 
two dimensional Ising system encountered here is to use ipi = (Pi Pi ) as the mean field. But once we treat 
the Ising spins in this way, nothing is really gained by introducing a source for the matter bilinears, and so 
we will only source the spin bilinear, repeating the mean field treatment of the previous section using this 
new order parameter. In the rest of this section we describe the results of doing this, and we compare them 
to those of the more traditional mean field treatment of the previous section. Although this second study 
will not be different in spirit from that of the previous section, it will give us some idea of the robustness of 
the mean field method. 

We introduce the source by inserting eS x t P fc P fc into the original path integral (15), and we define 
W(x) by calling the resulting path integral e w . Then our effective field is tf? k = dW/dxi and the effective 
action is 



A(<p) = W( X ) 



kj 



(45) 



As in Section 4 we define the mean field approximation by doing the spin sum with the spin dependent 
terms in the matter action fixed to their expectation values which are then self-consistently determined. In 
contrast to section 4, the mean fields here are the bilinears (P^ 1 ' P k ) and the derived expectations are now 
the linear and cubic ones (P k ) and (P[~ l PIpI +1 ) ■ 
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For the fishnet spin pattern p J k — tp,tp,tp' ,tp' ,tp,tp, . . ., there are three distinct values each for the linear 
and cubic cases. As explained in Appendix A. 2, where these expectations are evaluated, we denote the linear 
ones by </>, cj>' ', and <f> and the cubic ones by (+++), ( — h+), and (h — \). The matter integrals are evaluated in 
Appendix C.2, basically by substituting a new set of coefficients for the a, (3,^,6 coefficients of C.l. Now 
the effective string tension depends directly on the background fields 



T c s = — \ 



(At + ip + tp') 



a y 2(tp + tp') + tptp'/e)' 
As before the tension goes to zero as e — > unless tp' = 0(e). 



(46) 




■- l,e = 10-\p = 



y = 0.2 

V = 0.3 

iff = 



Figure 10: The spin energy per site for the massless case, p = 0, e = 10 -5 . On the left are the g 2 — 1 spin 
energy curves from Fig. 4 plotted versus (++). On the right are the analogous curves from the mean field 
treatment of Section 5. 

The energy function £(tp, tp') is certainly a different function of its two variables than the energy function 
of Section 4, both because the variables are differently defined and because the different structure of the 
source terms leads to an effectively different treatment of fluctuation effects in the corresponding mean field 
approximations. To display these differences graphically we choose a typical coupling g 2 = 1 for which they 
are evident. In Fig. 10 we show on the left the spin energy curves for this case from Fig. 4 in Section 4, 
replotted as a function of (++), the physical quantity measured by tp. On the right we show the corresponding 
curves from the treatment of this section as functions of tp. The most dramatic difference compared to Fig. 4 
is the move of the location of the minimum energy toward the origin. This effect is clearly accounted for 
by the change of independent variables from <j> to (++) . Next the shape of the curves on the left is slightly 
different in detail than those on the right, but the ordering of the curves is the same in both figures in 
the region near the minimum. Note that the <f> — <p' (tp = tp') curves in the left and right are essentially 
identical. However there is a significant lowering in energy of the constrained curves on the right compared 
to those on the left. This quantitative difference in the mean field treatments is reasonable considering that 
the constraints imposed on bilinears <p, p' in the present section are less severe than the constraints imposed 
on <j)' in Section 4. Hence the constrained system in the former case should be closer to the unconstrained 
ground state and the corresponding energies should be lower. Also the slicing of the two dimensional surfaces 
is slightly different in the two cases. This interpretation is supported by the fact that the ground state curve 
is not significantly altered. These two differences, the move toward the origin, and the lowering in energy of 
the constrained curves but not of the ground state curves are present in the spin energy for other couplings, 
but there are no other significant differences, so we don't display all those cases. 

Finally we turn our attention to the total energy per site. In Fig. 11 we display the total energy curves 
for the massless case and the same set of couplings used in Section 4. These graphs should be compared to 
Fig. 6. The most dramatic apparent difference is the movement of the minimum toward the origin as the 
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coupling increases. As discussed in the previous paragraph, this is explained by the change of independent 
variables from <f> to (++). For the strongest coupling displayed the minimum is for <p m .01. We also see the 
trend that the constrained curves generally lie lower in energy in Fig. 11 than in Fig. 6. But note that these 
constrained curves are not really comparable at strong coupling because in Fig. 11 most of the ip' choices are 
larger than the location of the minimum, in contrast to Fig. 6 where all of the 4>' choices are smaller than 
the location of the minimum. From Fig. 11, we see that the minimum energy is on the <p = <p' curve and the 
location of this minimum point depends upon the value of g 2 . For weak coupling, the minimum point tends 
to <p = 0.5 and it approaches ip = for strong coupling. The massive case p = 1 is presented in Fig. 12 and 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.05 0.1 0.15 0.2 0.25 0.3 

Figure 11: The total energy per site for the massless case, p — 0, e = 1CF 5 . 

should be compared to Fig. 7. In the latter the location of the minimum is near <j) = at weak coupling and 
approaches <p — 0.5 at strong coupling. In contrast in Fig. 12, again because of the different variables chosen, 
the minimum approaches <p = in the strong coupling limit in addition to the weak coupling limit. But 
the behavior of the constrained curves reflects the qualitatively different physics of the two limits. There is 
ferromagnetic behavior at weak coupling, with ip strongly correlated with ip' . At strong coupling there is very 
little if any such correlation. In summary, although the different variable choice has caused some dramatic 
differences in the energy curves between the two versions of the mean field approximation, the actual physics 
contained in the two approximation schemes is actually quite similar. And the basic conclusions we have 
drawn are robust at least against different styles of mean field approximation. 
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Figure 12: The total energy per site for the massive case, p = 1, e = 10 . 

6 Conclusion 

The mean field study undertaken in this article goes some way toward illuminating the role played by the 
fishnet diagrams in the sum of planar diagrams. For the <I> 3 theory studied, we identified a mean field 
description involving two mean fields <p, <fi' such that the fishnet diagram is singled out by the field values 1, 
respectively. But, within our mean field approximation, over the whole range of couplings the lowest energy 
is reached when <j> = <f>' and the effective string tension for these field values is 0. This indicates that the 
role of the fishnet diagram is negligible in the sum of planar diagrams of this theory. Since the underlying 
theory is unstable, this negative conclusion is probably reasonable. 

On a more positive note, our mean field method gives a concrete derivation of the kind of string description 
of field theory envisioned in the AdS/CFT duality of Maldacena. For generic values of the fields, we identify 
an effective string tension T e ^{(f), <j>') that characterizes the excitations close in energy to the energy associated 
with those fields. For <i> 3 theory it happens that T c ff = for the lowest energy configuration of fields. But in 
a theory like QCD which is supposed to confine quarks, we can easily imagine that the tension won't vanish. 
Application of our methods to this case will clearly be an interesting next step. 

But the well-known deficiencies of mean field theory still cast doubt on the reliability of such studies. We 
should regard them only as rough intuitive indicators of the physics potentially contained in these theories. It 
is therefore very important to develop calculational methods without such defects. One promising approach 
is Monte Carlo simulation. The worldshcet path integral for the q variables and the Ising spins clearly 
involves a positive definite integrand. The Grassmann be variables, introduced to give a local description, 
clearly introduce minus signs and ensuing complications for Monte Carlo methods. On the other hand, if one 
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integrates out the ghosts, positivity is again restored at the expense of some non-locality. Because the ghost 
dynamics is independent on each time slice, the relevant determinants are only one dimensional, and there 
is hope that their evaluation won't be prohibitively time consuming. Another numerical approach that is in 
principle exact is the block spin renormalization group. The feasibility of such numerical methods is under 
active investigation. We hope that the mean field results of the present article will be useful in interpreting 
any results that come from them. 

Even within the narrow application of mean field theory to the worldsheet for $ 3 theory, we have left 
some issues to be resolved. First, we left the alternative method [8] of approximating the matter variables 
with mean fields unexplored. In particular this method leads to a two-dimensional Ising spin system (41) 
which we despaired of solving analytically. However Monte Carlo methods are ideal for such systems and 
their application to this spin system is worth studying. The resulting hybrid of mean field methods and 
Monte Carlo methods provides a viewpoint midway between the intuitive mean field approach of the current 
article, and a more exact numerical approach. We have also proposed a worldsheet model for conventionally 
compactified <I> 3 theory in appendix D. The methods of the current article have yet to be applied to this 
model. The dimensional reduction of this compactified model gives a worldsheet model of two dimensional 
scalar field theory, which has already been investigated using discrete light-cone Hamiltonian dynamics 
(DLCQ) [16]. It will certainly be interesting to compare the results of this work with those of the path 
history worldsheet approach described here. 
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A Spin Sums in Mean Fields 
A.l Mean field is linear in P 

We turn first to the case where the mean field is taken linear in P. Then after replacing various monomials 
of P's in the action by their corresponding expectation values, we are left with the spin sums 




K J P' J 

K k^k 



1 - 44" 



In g 



s=±l 



(47) 



One option is to also replace the s k = 2P k — 1 by 2<j) k — 1 as discussed in the text. But when a simple pattern 
is assumed for the K k , the sum is just that for the partition function of a one dimensional Ising model with 
uniform magnetic field and can be done exactly. 

k is uniform, the spin sum is just the Nth power of the 2x2 matrix 



When K k 



1 



whose eigenvalues can be immediately written down 
1 + e K ± - e K ) 2 + Ag 2 e K 



t± = 



= e K/2 |^cosh 2 ± Y sinh2 2 + ^ 



(48) 



(49) 



Thus there are two branches for the spin contribution to the energy per site — hxt±(t£), the branch with t + 
clearly being the lowest energy. Altogether the k dependent terms in the exponent of the Boltzmann factor 
are MN[lnt + — </>k] so the saddle point in k determines its relation to <p 



d_ 

dn 



lni+ 



sinh(«/2) 
^g 2 + sinh 2 («;/2) 



(50) 
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2^(1 " 0) 1 2^(1-^) V ^(W) 



(51) 



The second equation enables expressing the spin energy per site as a function of <p, 

£ s = n(j)-\nt + (52) 

^ J \ 2^(1 V 4 <KW) J 1 2^(1 -<^) V ^(1-^) J 

For the fishnet spin pattern, which repeats only after four time steps, we have a more complicated matrix 
to diagonalize. To force the 0', </>, <j>, <f> pattern for the <^'s, we need the K J k in the pattern k', k, K, k so the 
transfer matrix is a power of the 4 step matrix 

e K ge K \ / e K ge K \ ( e K ge K \ ( e K ge K 



T *-U i JO iJU i J- ,54) 

The eigenvalues of a general 2x2 matrix can be written 

4 = TrT± ^ T f 2 - 4dc ^. (55) 

It is straightforward to work out the trace and determinant: 

detT = e 2K+K+K ' (g 2 - l) 4 (56) 
TrT = e 2K+K+K ' + g 2 e 2K+K ' + 2g 2 e K+K+K ' + g 2 e 2K+K + g 4 (e 2K + e K+K ') 

+2g 2 {e K+K + e K+K ') + 2g 2 e K + g 2 e K + g 2 e K ' + 1. (57) 

After a bit of algebraic manipulation the eigenvalues can be written 
1 



(! + e (2«+K+ K ')/2)2 + 2 52 e «(l _ e (K+ K ')/2)2 + ^4( eK + p (K+ K ')/2)2 + - g 2^K + + e «) 



2 



±J (1 - e (2«+^+«')/2)2 + 2g 2 e K (l + e (^+«')/2)2 + g4( e « _ e (K+ K ')/2)2 + ^2( e K + e n')(l + e «)2 



2 



(58) 



Again it is clear that t\ is the branch with lower energy. The exponent of the Boltzmann factor is now 
MN{\nt\ - 4>{K + 2k) - fin') /4, so we have 



v / (TrT)2-4detT 2 ^(TrT) 2 - 4det T 
e K ' (e 2K+if + ,g 2 e 2K + 2g 2 e K+K + g 4 e K + 2g 2 e K + g 2 ) - tt 



sj (TrT) 2 - 4 detT 

1 

2' 



(59) 



(P) = -9 K lni 4 .=^lni 4 _ (60) 



e K ^ e 2 K+K ' + ~2 e 2 K + 2g 2 e K+K ' + ,g 4 e K ' + 2.g 2 e K + .g 2 ) - t 4 
x : 1 1 / l.uTT ~ 
The equality of the two expressions for cf) determines K: 



(61) 



e K = eK e^+9y + l + e^ 
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We also list the expectations of the higher monomials in P that occur in the action. Let (k\, K2,k^, K4) 
denote (k, K, k, k') in a general cyclic order. Then, in this notation 



TrT 



d Kl TvT 



„Kl+K 2 +K3+K4 



+ 1 + g 2 [e* 1 



+K2+K3 



e 



K2+K3+K4 



e «2 + K3 _|_ gK3 + «4 _|_ g/<2 



+e 



e Kl+K2+K3 + «4 _|_ ^2 j~gKl+K2 + K 3 _|_ g 



«4l _|_ g 4 ( e K i+ K f 



K1+K2 

K2 + K4 



gKl+K 4 _|_ gKi 



e 



K1+K3+K4 



_|_ gKl+K2 + K 4 _|_ gKl+K2 



+ ^4gK 1+K3 



1 1 ; 1 1 



= § [ e «2+«3 + «4 + g«3+«4 + g«4 + ]_ + ^gK 2+K4 + g« 2 + « 3 + g«2 + g«3)] _ 

Then we obtain for P's on consecutive temporal sites: 
1 d Kl TrT — TrT/2 



(P 1 ) 
(P 4 P 2 ) 



2 ^(TrT) 2 -4dctT 

„Kl+K 2 



(63) 
(64) 

(65) 
(66) 



i 4 
5 + 



+ .g 2 (e K3 +e K4 + l))(P 1 ) 



(g«3 + «4 + ^2(g K3 + g«4 + ^ + 



g(e K 



3+K4 



+ g 2 e K * + e K3 +l)Y 



V(TrT) 2 -4dctT 



v /(TrT) 2 -4dctT 



(67) 



Ni = 



-- ( e K 3+^4 + ~2(gK 3 + gK4 + ^ + g/Cl+ZCa + Ks+K*^ + g«3 + «4) 



+.g 2 [e Kl+K2 (l 
+ ( 1 + e «3+«4 . 



= K3 + K 4 



)(1 



= K3 + K4 



') + e' 
) 



K1+K2 + K3 + K4 /„K3 



(e K3 +e K4 ) 



_|_gK3+K4^gKl+K4 _|_ gK2 + K3^ _|_ gK 3 + K4^ -|- e K3 + K4 )(e Kl + e" 2 )] 

+.g 4 [e Kl+K2 (e K3 + e K4 )(l + e K3 + e Ki ) + e Kl+2K3+K4 

+ (1 + e K3 + e K 4)( e K l+ K 3+ K 4 -)- g K l+ K 4 _|_ gKl~j _|_ g«4^ -|- g«2 + K3 + K4 _|_ gK4 _|_ gK3 + K4^ 

+ (1 + e K3+K4 + e K3 )(e K2+K4 + e K2+K3 + e K2 + e K3 )] 
+g 6 [(e Kl+K3 + e K2+K4 )(l + e K3 + e K4 ) + e K3+K4 ] 



(P 1 P 2 P 3 ) 



gKl+K2+K3 
r + 

gK!+K2+K3 



(.9 2 + e K4 ) 



(3 2 



(P 1 ) 



g (e K4 + 1) y 
v/(TrT) 2 -4dctT 

v/(TrT) 2 -4dctT 



(68) 



(69) 



N 2 = 



^2 _|_ p K 4 j ^7" _|_ gKl+K2 + K3 + K4 _|_ ^6gKl+K 3 + K4 



+g 2 |^2e Kl+K2+K3+K4 + e Kl+K3+K4 + (e Kl + e K3 )(p K2+K4 + 6 k 2+2k4 _|_ e 2«4 _|_ e K4 )] 

+g 4 [e K2 (l + e K4 ) 2 + (1 + e K4 )(l + e K2 )(e Kl + e K3 ) + e Kl+K3 (e K2 + 2e K4 )] (70) 

By inspecting Ni and -/V2, we see that the double P expectations (P l P 2 ) are symmetric under (ki,«3 



<— > 

1 p2\ 



(«2, K4). This means that the fishnet patterns + + H — and + H h give identical values, (++), of (P L P 

and similarly the patterns H 1 — h and — H — h give identical values ( — h). The triple P expectations have the 

symmetry n x <-> k 3 which implies identical values, ( — h+), for — +++ and +H h But the triple expectations 

for H h+ and + + H — arc unrelated d, and we denote them by (h — h)and (+++) respectively. Thus 

for the fishnet spin pattern (P 1 P 2 ) assumes two distinct values and (P 1 P 2 P 3 ) assumes three distinct values. 

For k! = k, it immediately follows that K = k and <\>' = <p. Therefore this case should reduce to the 
uniform mean field results. Indeed, one finds t\ = (t+) 4 as expected. Then it is straightforward to show 
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that the expression for <f) reduces to (50), and the expressions for the higher correlators become 



(p'p 2 ) £1-*+^! = ^! (71) 
t-y. t-\- — t— t-\. 



(P^P 3 ) - " i±_L=^(|_) . (72) 



i 2 + i+-t_ 

As discussed in the text, the e — > limit is smooth if one of the mean fields also goes to zero. For the 
fishnet pattern we achieve this by putting <fi' = ef and taking the limit at fixed /, <f>. We therefore need to 
specialize the spin contribution to the energy per site to this case. We note that <f>' —* requires k' — ► — oo. 
Holding cj) away from its endpoints 0,1, we find in this limit 

-+ £gO±g) (74) 

4 -» TrT = (1 + .gV) [1 + 2.g 2 e K + ,g 2 e 2K ] (75) 

<++> - <"> 

<+++> - (7^ (77) 

, 2 K > 1 + 2(1 + g 2 )e K + (2 + (1 + g 2 ) 2 )e 2 " + 2(1 + g 2 )e 3K + e 4 * 

* 5 6 (1 + 2.g 2 e- + 5 2 e 2 «) 2 1 ' 

/ x k'j. l + e K + 9 2 e K + e 2K 

( - +) - ^TTWW (79) 

K , <Pe K l + e K +g 2 e K + e 2K 
\ ++) -» e 1 + gK 1 + 2.g 2 e« + .g 2 e 2 « ( j 

, 0e"(l + e") 

(+-+> ^ e TTW# ( } 

Notice that the quantity ( — h), which figures in the effective string tension is proportional to <j>' in this limit: 
/ +) _ e^l + e^d + e^+g^ + e^) = 

^ 1 M + 2(l + g 2 )e K + (2 + (1 + ,g 2 ) 2 )e 2 « + 2(1 + g 2 )e 3 « + e 4ft ~ v 1 j 

For visualization purposes we display graphs of the quantities (++), (+++), and G in Fig. 13. 

Moreover, we see that in this case e K can be eliminated in favor of <f> by solving a quadratic equation: 



e 



V(l -20) 2 +40(1 -0)/.g 2 - 1 + 20 



2(1™0) 
Then the spin energy per site, 



(83) 



£ s = 1 [0V + 2n<f) + K<j)-\n t%] (84) 

- J [(2k + X)0 - ln(l + $V) - ln[l + $V + $V(1 + e K )]] 

- i[(l-0)ln(l-0) + 01n0-01n.g 2 ] 



can be expressed as an explicit function of 0. Note that the term 0'k' = O(elne) and therefore drops out 
for e -> 0. 
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A. 2 Mean field is quadratic in P 



We now turn our attention to the case where the mean field is taken quadratic in P. Like the previous 
subsection, we replace the P, PP's in the action by corresponding <f) and ip, then we are left with the spin 
sums 



E ex p i E 

k,j . 



(86) 



For the uniform case in which xl = X: the spin sum as before is just matrix multiplication of the N 2 x 2 
matrix 



e x g 



which has eigenvalues 



t± 



9 1, 

1 + ex ± -e*)2 +4g 2 



(87) 



(88) 
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For large N, t+ is the largest eigenvalue of the transfer matrix. For our purpose, we need to calculate 
expectation values of single, double and triple P. After some algebraic manipulation, one gets 

4> = (P) = <*-W(l-e*g + f£ (89) 
2^(l-ex)2 + 4g 2 

v = (pp<) = ^-i + V(T^VTW) 

(1 + ex + - e x)2 + 4g 2 )^(l - e x) 2 + 4g 2 
(PP'P") = 2e 2 MeX-l + N /(l -ex) 2 +4g 2 ) 

(1 + ex + ^(1 - e x) 2 + 4.9 2 ) V(l - cX ) 2 + 4 .9 2 
For the fishnet pattern we have to evaluate the spin sums by diagonalizing the 4 time-step transfer matrix 

f OCT ?)(* 0(f 0- 

To present the results of this analysis in a reasonably compact way it is useful to rewrite T4 using Pauli 
matrices, i.e. 

and a similar definition for x'i- Then 

T ±=(^~) {I + U ' ■ (T )( I + U - (T )( I + U - (T )( I + U ' ■*) ( 94 ) 

One can explicitly calculate u and u' 

( 2g e*-l\ , ( 2g e x ' - 1\ 

«= h-^r. 0, — — , u'= — ^-7, 0, 1 



1 + ex ex + 1 y I 1 + ex ex' + 1 / 

^,(1-^+5, (95) 



(l + ex) 2 ' (1 + ex') 2 

It is straightforward to show that the eigenvalues of T4 are 

s 2 / 2 

l + e*\ / 1 + e 



4 = (j-Y-) \ ^~) (\ u ' + u \ ± \ /l + u2u ' 2 + 2u - u ') ( 96 ) 

= \ (^ v /.9 2 (2 + ex + e x')2 + (i _ e x'+x f 

±^jg 2 (ex + e x')2 + 4 ^2(! + ^2) + (1 + e x+x')2"j (97) 



To simplify our formulas we introduce some notations: 

A=u + u', B=l + u 2 u' 2 + 2u-u', Q = (1 -u' 2 )A+ A 2 u', g = (1 -u 2 )A+A 2 u (98) 
With these notations, we can rewrite the transfer matrix as follows 

Ta -\-^T) I — 2 — J I 2& A 2 +B-2G 3 (99) 
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and the eigenvalues can be written as 



2 



2 



t± = {^f) [^j {j? + B±2^S) (100) 

Notice that A., B and t± are symmetric under \ ~ * x' while Q — > (?. However, one can easily show that 
Q ■ Q = Q ■ Q = j4?B is symmetric under \ — * X 1 ■ F° r our problem, we need to calculate correlation function 
for (P), (PP') and (PP'P"). The results are 

, g 3 + Vg~g g 3 + Vq~g ~_ g 3 + VG~U ggi(c* -i)-2 g 2 g 3 

= eH^-Gz-Vgrg) g^grfl (102) 



8e*+*'[(g 2 + cx+x')(g 3 + + ggi(l + c*')] 

(^ 2 +fi + 2Ve _ a)Ve _ a(l + cx)2(l + cx')2 1 ; 

8c 2 X[(g 2 + c 2 X ')^g-g _ (g2 _ c 2 x) g 3 + 2 gg ie x'] 

(i + c x)2(i + cx')2(^ 2 + b + 2Ve~g)Ve~G 

8e 2 x' [(g 2 + c 2 x) - (ff 2 - c 2x )g 3 + 2 gg lC x] 

(i + ex) 2 (i + e x')2(^ 2 + B + 2Vg~g)Ve~G 



For large M, N the spin energy per site is 



£ S (X,X') = - J(ln4( X ,x') - 2^x - + ^ ln(l + + 0' + 20) (106) 



One can easily show that when X = x' (uniform case) 



ex-l + v /( e x-l)2 + 4g2 
2 v /( e x-l)2 + 4 g2 



¥> = 



e x - 1 - 2g 2 + ^/(ex - l) 2 + 4g 2 



2(e*-g 2 )V(e x -l) 2 + 4g 2 



r+ - t_ 



l + (2 + eX) 2 g2j [( 1 + 2 3 2)2 + 3 2 e 2 x]j 

-l + (-3+(l + e*) 2 ) g 2 + t+-t_ 
2(t+-t_) 



(107) 
(108) 



<- ++ ) - ^^-^V 7 ^^^) (109) 

(1 + ex + - e x)2 + 4g 2 ) 2 \/(l - e*) 2 + 4g 2 

(+-+> = (+++> = (-++> (HO) 
as expected in Eqs. (89), (90) and (91). In this special case, the spin energy can be written as 

l + e * + v /(l- e x)2 +4 g 2 d 
f s = (^X-ln S£ — ^ + -01n(l + p). ( 1U ) 

Another special case is ip' = or x' — » — oo. In this case 

t± - J (\/g 2 (2 + ex)2 + 1 ± ^(1 + 2g 2 ) 2 + g 2 e 2 x) 2 , (112) 



(113) 
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-l-(2 + e* (2 + e x )) g 2 + t+-t 

2(i+-i_) 
-1 + (-2 + e 2 *) g 2 + t+ - t_ 

2(t+-t_) 
e x <p [i+ (e x + 1) + .g 2 - e*] 
*+(*+-*-) 

In this limit, (-++) = (h — h) = and (+++) is 

e 2x 9 2 [g 2 + (l + ex) 2 g* + t + -t- 



(+++> 



(t+ - t_) [1 + 2.g4(t + - t _) + g2 (eX + 2)2(<+ _ 



In this case, the spin energy is 



v /.9 2 (2 + e x)2 + i + v /(i + 2 ff 2 ) 2 +ff 2 e 



2x 



+ \<px + ^(4> + <t>' + z<t>) Mi + p) 



(114) 
(115) 
(116) 

(117) 

(118) 



B The coefficients 

B.i 4 = {Pi) 

The Gaussian integrals over ghost and q fields for the fishnet spin pattern requires the extraction of the 
coefficients of the various bilinear forms. Let us write the ghost action for a single 6, c pair on time slice j as 

S g = a j ^b k c k + f3 3 Y^(-) kb kCk+l 3 ^{bk+i -h)(c k+1 -c fc ) 

k k k 

+S j ^HHbk+i - b k )(c k+1 - Cfc ). (119) 

k 

Then for four consecutive time slices (and periodic boundary conditions) we read off: 



a = 



P 1 = -P 3 = 



7 

5 1 = 

a 2 - 



- + 1 + p(l - 0') - 277(1 - 0') - r?0j (++) - (+++) - (+-+) + 2(0 + 0') - 
+ Q + 1 + p(l - 0) - 2r?(l - 0) - r?0') (-+) + 2/7(1 - 0)<+-+) 
+2r7(l - 0')(+++) + »7W + <f>'){-++) 
(-|-1 + P(1- </>') - 2r7(l - 0') - (++) + (+++> - (+-+) + 2(0 - 0') 
+ (j + 1 - P(l - 4>) + 2/7(1 - 4>) + (-+) - 2/7(1 - 0)<+-+) 
+ 2r7(l-0')<+++)+r?(0-0')<-++) 

= 7 3 = (i - 0)(i - 0') - - 0')«-+) - (-++)) - ^(i - 0)«++> - (-++)) 

= -<5 3 = 77(1 - 0')<-+) - 77(1 - 0)<++) - (0 - 0')<-++) 



(120) 



a 4 - 



((++> + (-+» - + l + p(l-0) - r?(2 - 0) +40(1-0) 



-2(-++) + 477(1 - 0)(-++) + #((+++> + (+-+)) 



(121) 

(122) 
(123) 



(124) 
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/3 2 

7 2 
S 2 



((++> - (-+)) ( -- + 1 + p(l - <f>) + r?(2 - 0) ) - #«+++) - (+-+)) 



7 4 - (1 - - ^(1 - 0K++) + (-+> - (+++> - (+-+)] 



S 4 = 7 7 «++) -<-+)- <+++) + <+-+» 



Similarly we write the bilinear form for the q's as 

& k 

with the ± sign alternating with time slice according to the pattern + H h H •••• Then 



2e 



2e 



(125) 

(126) 

(127) 
(128) 



(129) 



(130) 



B.2 tf = (W 



Just as the previous subsection we need to extract coefficients of the various bilinear forms. The ghost action 
and q action have the same form as the previous subsection. After reading off, coefficients are 



a 1 



7 1 



a = 



5 3 = 



1 



+ i + P (W) 



if- 



1 



+ 1+P(W) 



■'-<+++} -<+-+>} 



- + 1- P (W) 



<p' + 2(0 + 0') 

<p'-2(4>-4>') 



-<+++} + (+-+» 

7 3 = (1 -</>)( W), ^ = -^ = 



- + i + p(i-</>; 



2 



e 



- + l + p(l-0) 



¥?' + 40(1 - 0) - 2(-++) 



7 = 



7 4 -(l-^) 2 , £ 2 = -«5 4 = 0, 
for the ghosts and for the g's they are 



a = — , p ■ 



2e ' ' 2e 

Note that we denote coefficients correspond to (PP') case with tilde. 



(131) 



(132) 



C Fishnet Determinants 

C.l Single P case 

The easiest way to do the q integrals is to diagonalize the bilinear form in the exponent by going to normal 
modes. We describe this procedure by first going to normal modes in i, and then in j. We let the i labels 
range from i = 0, 1, . . . , M and choose Dirichlet conditions at the boundaries: q 3 — q 3 M = 0. Thus our system 
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carries zero transverse momentum (we are working in the center of mass system). Then the mode expansion 
reads: 

1=1 

For completeness, we also mention the mode expansion for periodic boundary conditions (which will be 
useful in our discussion of the ghost determinant.) 

-M-l 

1 V I r / /,■ ' 

(134) 



q k = V T7 E A l ex P { ^TT } ■ Pcriodic bx - 



Defining the frequencies uuf = 4 sin 2 lir/2M, we find 



E(^+i-^) 2 =E^V [£^m- ; <4, Periodic b.c.) (135) 
k i \ i / 

For the remaining terms in the exponent the g^'s have a different coefficient that the g^ en 's. Thus in 
addition to the usual orthonormality condition 



fc=l fc=o k ' 

we also need the corresponding formulas when the k sum is restricted to even and odd values. A short 
calculation shows that 



2 M_1 irl'k nlk 1 M 1 (2ni(l' + l)k^ 

MEH k Bin irS in ¥ = -^M- ! , -^R^P ' M = &l>M/2-l mod M, (137) 

fc=l fc=0 k J 

where the periodic case assumes that M is even. From these we derive the identities 

2 ^ . wl'k . irlk 1 , 

M L sm ^ sm m = 2^V + ^,m-0 (138) 

fe= odd 

i E -P{^f^) - (139) 

k=odd L - - J 

2 ^ 7T?'fc . 7r/fc 1 

M E Sln lw- Sln M" = 2 (J, '' , " 5, '- M - ,) (W0) 

ij E ex p — H7 — = r- (141) 



M ^ r I M 

/eleven 



The respective sums are then 

E = I E Q% + i E QtoL-v E «i a = J E Q? - \ E , (m 



k=odd I I k=even I 



for Dirichlet boundary conditions and 

E ^ = o E A i A M-i - ? E A ; a m/2-;; E = o E A \ Aj M-i + 9 E A i A M/2-i ( i43 ) 



k=odd I I k=even 
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for periodic boundary conditions. Depending on the value of j, cither (++) multiplies the even sum and ( — h) 
multiplies the odd sum, or vice versa. Recalling that a = ((++) + (-+))/2e and (3 = ((++) — ( — h))/2e,we 
find that the entire exponent for fixed j for the q integral can be written as —a /2m times 



(144) 



for Dirichlet conditions and 

^lA^M-x + "£(4 - M'^m-i Ai M-i) ± /?£(4 ^'" 1 )(</ 2 - i - ^m/2-;)- (145) 



for periodic conditions. The ± stays fixed for two time steps and then switches. For example the simple 

function \f2 cos[7r(2j — l)/4] produces the signs H h H 1 as j = 1, 2, 3, 4, • • •, so we can write the 

complete exponent as ~a/2m times: 



5>?Qf + a£(Q? QF 1 ) 2 + /?£ ^cos J(2j - 1){Q\ Qf^u-i - Q'm-i)- 



(146) 



1,3 1,3 1,3 

with a similar expression for the periodic case. 

Now we are in a position to choose normal modes in j. Let us let the j values range from to N, and 
impose Dirichlet conditions at j — 0, N: Q° = = 0, so that wc have for the Dirichlet case 



2 



(147) 



n=l 



Defining the frequencies uj'^ = 4 sin 2 mr/2N, we easily find 

j2Qf = E^r 2 ; - of- 1 ) 2 = £^e (W8) 

j n j n 

But we also need the more complicated sums in the last term of the exponent. To do that we need 



2 

N 



£ cos ? j-. 



i=i 



. mr? . n7r(i — 1) 

AT 



. n'irj . nir(j — 1) 
sin — - - sin 



N 



N 



(149) 



, . nir , n'lr 2 -r-^ n ( 1\ mr / 1\ n'n ( 1\ 
4 sin — — sin — — — > cos — 7 — - cos — ? — - cos — — ? — - . 150) 
2N 2N N 2 \ 2) N \ 2 J N \ 2 J K J 

We shall restrict N to even values, in which case we find that this expression is essentially block diagonal 
with the blocks being at most 4x4: 



mr n'n 2 7r / 1 

4 sin — — sin — — — > cos — I j — — 
2N 2N N ^ 2 V 2 

.7 = 1 V 



COS 



mr 



1 



COS 



n it 
~N~ 



1 



„ . mr . n 7r 
2 sin — - sin — — 

27V 2JV 



S n ,N _ n + S n , N +n + 5 n , n _N - 5 n , 3N_ 



(151) 



The exponent now becomes —a /2m times 

+ a^)Af + /?^£ A l A M-l"Wn' k',f -„ + <kf +n + K>,n- 



n' n 



£(w? + +/3V2 E K',f + <V,f +n + <Wf ~ <W 

2,ra l<M/2,n 



(152) 
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For I < M/2 the mode A? couples to itself as well as A N J 2 7 n , A N J 2 + n for n < N/2 and to A^Z", ^m_?~" 
for n > N/2. For clarity, give a different letter name to each of the four possibilities: 



( A? for I < M/2, n < N/2 

Bf~ n for I < M/2, n > N/2 

CM-i n for I > M/2, n < N/2 

. D l7 N / 2 for 1 > M A n > N/2 



(153) 



Then the bilinear form can be written 



KM/2,n<N/2 

+[3V2 [A?«,' n {C?u/ N/2 _ n + Dfu' N/2+n ) + B?w' N _ n {C?u,' N/2 _ n D?u/ N/2+n ) 

KM/2.n<N/2 



(154) 



The Gaussian integral then involves the determinant of a block diagonal matrix, the size of each block at 
most 4x4. The computation of each determinant is not hard: Each 4x4 block has the structure and 
determinant 



M 



a ef eh 

b gf -gh 

ef gf c 

, eh —gh d 



(155) 



detM = abcd + 4(efgh) 2 -ac(gh) 2 -ad(gf) 2 - bd(eff - bc{ehf . (156) 
For the 4x4 matrix for mode (I, n) with I < M/2 and n < N/2, we have the following values for a,b, . . . ,h: 



a = u)f +awf; 



b = uf +awj?_ n ; c = ij? m _ x + a^ /2 _„; d = uj 2 m + aJ 2 



ef 



N/2-n _ 



V2 ' V2 
Plugging these values into det M we get 



N/2+n . P U1 N-n U> N/2-n 

; gf = - 



V2 



gh 



M-l ' lxuJ N/2+n 

71 



(157) 
(158) 



dctM ; ,„ = 16(a 2 -/3 2 ) 2 sin 2 — cos 2 — + 16aK-/3 2 )K 2 cos 2 — +^_ ; sin 2 — ] 



+4a 2 (, 



n7r 



cos" — + lOm_! sin" ) + 64 ( a 2 — 



N 



N J 



ft 2 



+ a sin 



M 



16 sin z 



Ztt 
M 



In obtaining this result we took advantage of the explicit forms for the id's: 



u t = 

, ,' 2 

%/2-n — 

In particular note that 



, ■ i 2 
4 sin — , u M _ 

4 sin 2 (1-™) 
V4 2NJ 



2 

4C ° S 2M' - 



12 



nir 



4sin^— , ^_„ = 4cos^ 2iv 



nir 



"N/2+n -4COS- 



/7T 7J7T 

V4 ~ 27V 



2 , 2 



UJ' 2 + LU' 2 



u?u 2 m-i = 4 sin 2 



In 

M' 



UJ N/2~n + ^W/2+n ~ 4 
12 12 t • 2 ^n /2/2 t 2 ^n 

w„W;v- n = 4sin — , u N/2 _ n u N/2+n = 4 cos 



N 



N 



(159) 

(160) 
(161) 

(162) 
(163) 



Note that detM; jn is a quadratic polynomial in sin 2 (n7r/iV) which we can factorize by finding the roots r±: 
detM,, n = 16( a 2 -/? 2 ) 2 (r + -sin 2 ^)(sin 2 ^-r_). (164) 
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We can conveniently parametrize the roots as 
r + = Ri cosh 2 m , 

Ri 



r_ = — Ri sinh k; 



l a 3 - a[3 2 + a 2 lir 

1 + 4 W^W? cos m 



Ri sinh m 



l Rf + "VA-i + 2a)[(u;? + 4a)(g;^_ t + 2a) - 8/^ _ , 



4(a 2 - /3 2 ) 2 



The result of doing the integral over all the q's is then simply 

'2nm\ dMN/2 



a J 



]} det~ d / 2 M hn . 

KM/2.n<N/2 



(165) 
(166) 

(167) 



(168) 



Note that for M odd, which we can take for convenience, there is no special unpaired I = M/2 mode. 
However, since we must restrict N to be even there is a special n — N/2 mode which must be handled 
separately. However, as far as the energy is concerned, it's presence or absence doesn't matter because the 
energy is extracted from the coefficient of N in W as N — ► oo. 

Using the above factorization, we can perform the product over n by using one of the identities: 



N-l 



2 n7T , . ,2, 

hiii — + 4 sinh £ 



n( 4 si 

n=l 
V-l 

J] (4cosh 2 e-4sin 2 ^) 

n=l 
V-l 

J] (4cosh 2 e-4sin 2 ^) 



4 cosh' 



c n (^i 



n<N/2 



-in" — + 4 sinh 



(169) 



sinh z £ 

4 sinh 2 £ TT f4cosh 2 e-4sin 2 ^ N ) 2 = ^^, for N cver<170) 
■ LX V W / cosh £ 



n<N/2 



TT f4cosh 2 ^-4sin 2 ^) 2 = ^L^ j for7Vodd (m) 
V jv / cosh £ 



KN/2 



Since we want to product over only the modes with n < iV/2, we will require the square roots of these 
identities: 



T [ ^4sin z f 4 siulr 



n<N/2 



N 







JJ (4cosh 2 C-4sin 2 ^) 



sinh AT£ 
sinh 2£ 

sinh 7V£ 
sinh 2^ ' 



n<JV/2 

which are now specialized to even N. the appropriate £'s to use are given by 

sinh £}_ = \J —r- = \f~Ri sinh m , cosh = y 7 ?^ = c °sh K z > 
which can be solved: 



£L = In ^v/^sinliK; + \j I + Ri sinh 2 K^j 
£ l + = In (^s/R~i cosh/i; + ^/ Ri cosh 2 — 1^ . 



Thus we have 

(2fm/ fl )- JM/2 e ff ' 



n n detM,, n = n 

KM/2n<N/2 KM/2 



2 JV _ 2 sinhiV^ sinhiV^ 



-, -d/2 



sinh2^ sinh2£< 



(172) 
(173) 

(174) 

(175) 
(176) 

(177) 
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From this we extract the q contribution to the energy per site (£ q = —W q /MN): 

d 



2M 



E &+d+Ma 2 -/3 2 )]. 



(178) 



KM/2 



We will eventually need the contribution of the general ghost action which will have the coefficients 
of the terms (b k+ \ — b k )(c k+ i — c k ) different for k even and odd. Unfortunately for Dirichlet boundary 
conditions this translates to non-local expressions in mode space. However, the bulk energy, which we have 
been extracting is the same for Dirichlet and periodic boundary conditions, and for even M this more general 
anti-ferromagnetic pattern of couplings is local in mode space: 

a^2b k c k + (3^2(-) k b k c k + -f^2(b k+1 - b k )(c k+1 - c k ) +5^(-) k (b k+1 - b k )(c k+1 - c k ) = 



aB C + (a + k))B M/2 C M /2 + ^ B[C M -i 

i^0,M/2 

M/2-1 . . M/2-1 



a + 47 sin — 
M 



si 

1=1 ^ 



In lir 

sm TT cos TT 

M M 



+ 0(B o C M/2 + B M/2 C o ) (179) 

E b m/2+iC'm-i ((3 + US sin cos ^ 
1=1 ^ 



The corresponding ghost determinant for a single time slice is easily read off in this periodic case: 

M/2-1 r 



Detf - {a 2 +A ai -0 2 ) J] 

1=1 

= (a 2 +Aa 1 ~0 2 ){ 1 2 -8 2 ) M l 2 
-- In 



a 2 -0 2 + Aaj + 16( 7 2 - S 2 ) sin 2 — cos 2 ft 
v 'MM 



/ a 2 - /? 2 + 4a7 



4( 7 2 - <5 2 ) 



+ 4 I' 



.! /sinhMCp/2 
\ sinh Cp 

a 2 - 2 + Aaj 
4( 7 2 - <5 2 ) 



(180) 
(181) 
(182) 



The prefactor (a 2 + Aoff — 2 ) shows that the determinant is zero for 2 = a 2 + Aaj which is due to a zero 
mode that occurs in that case. But this is insignificant for the bulk energy per site: 




S 2 + 



a 2 ~0 2 



+ «i7i) ="2ArE& ( 183 ) 



Here j labels the time slice. For the fishnet spin pattern, there are two distinct types of time slice and two 
different values for , which we simply call £1 and £2- Then we write the total matter contribution (q plus 
ghosts) to the energy per site, including the boundary term for the massive case, as 



£m(0,0') = 



d3<f>- 



■mi+p)-^i+6)+^ E ie + +e-+M* 2 -0 2 )}- 



(184) 



KM/2 



In view of the complexity of this formula, it is useful to explicitly display simplifying special cases. First 
take the uniform case 4>' = 4>. Then = J g — 5 J g — 0, and 



(++) i +/HJ -o) - „(2 -<>) ) 4 2o(1 -<>) - i - )(1 - „(2 -o)) 

(++> 
e 

7 =(1 



7,(1 - 0)«++> -(+++», 



<++> 



(185) 
(186) 
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Then the matter energy per site simplifies to 

M-l 

2-r-v- ■ r, • M 



^ ln (l + P) + 4 E ln 
1 



e . 2 In 

(++y sm 2M 



"2 In 



«7e a z 



(++> 



T + tT 



eery a 2 



(++> 2 (++> 4 



^ ln(l + p) + — Re <JiLi 2 

2 7T 



(++) 



^ln 
2 



£21 4 + L 



eoey a 2 



++ 



(++> 



(187) 



The second special case is one in which we take e — > holding / = <///e fixed. In this limit a 2 — /3 2 



(++)(-+)/ ( 2 = fG(++)/e is of 0(l/e) and so are the a 2 — 0\ . Thus the £ ff simplify to 



(188) 



Keeping only the leading terms as e — > wc find 



In- 



+ ln y/2 + fG + (++)( P - v) - - (+++» 



+ ln + JG + (1 - 0)(2 + ( ++ )(p - v )+ v (+ + +)) (189) 

The simplification of the q determinant is less dramatic. We list the limiting forms of the various ingredients: 



Ri 



Ri sinh Ki 



1 + 2/G Itt 
1 + -PG^ C0S M 



'i? 2 + sin 2 ln 8(l + /G)-4sin 2 Z7r/2M 



2M 



J2 G 2 



Ri 



ln(a 2 - /3 2 ) = ln 



fG(++) 



(190) 
(191) 
(192) 



Combining these ingredients we end up with total energy per site for <f>' = ef: 

d 



f,(fl = ^ln(l + P) 



2M 



2 K/G) + ^+d] 



KM/2 



In V2 + /G + (++)(p - t?) - 2r ? ((++) - (+++)) 

In V / (++> + /G + (1 - + (++>(/> - i?) + v(+++)) 



(193) 



We note that the / — > limit of the right side is finite because ln(/G) + £ l + — ► ln(2^/cos(7r//M)) in that 
limit. 



KM/2 



sin 



2M 



1 + cos 2 



2M 



' cos — + sin — — - 

M 2M 



1 + cos 2 



Itt 
2M 



+ - A ln 2 - d - ln v/2 + (++> (p - t,) - 2t 7 (( ++ ) - (+++» 
-| In v/(++) + (1 - m + (++)(/> - ri) + ??(+++)) 



(194) 
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C.2 Double P case 



The above computations are general so one can obtain the results for the double P correlator by replacing 
a — > 5, etc. Instead of repeating the same lengthy procedure as before we just quote some important results. 
The q contribution to the energy per site is 

d 



where 



£q 2.1/ 



]T [? + +d+ln(5 2 -/3 2 )], 

KM/2 



£_ = In yyRi sinh Ki + y 1 + Ri sinh k; 
£ l + = In (^fRi cosh«; + \j Ri cosh 2 m — l^j 

Ri 



„ a 3 - a/3 2 + a 2 Ztt 
1 + 4 cos ■ 



(52 _ £2)2 M 



i?; sinh K; 



Ri cosh = 



Rf + 



4(5 2 -/3 2 ) 2 



ft 



u^-t + 25) [(^ 2 + 45)(^-z + 25) - 8/3 2 
4(5 2 -/3 2 ) 2 



The ghost energy per site is 



27V ^ 



'5 2 -/? 2 



5?"/3 2 



+ a i7i 



~(a+6) 



(195) 

(196) 
(197) 
(198) 

(199) 
(200) 



(201) 



here we have use the fishnet spin pattern where there two only two distinct types of time slice and two 
different values for and we will call them £i and £2- The total matter contribution to the energy per site 
as 



M^') = -^i+6) + ^ E [^+d+ln(5 2 -/? 2 )] 



(202) 



KM/2 



There are two special cases. First we take the uniform case in which ip = ip'. In this case, J g = 0, (3 = 0, 
5 = ip/c and 



op 



<p- (+ + +> +2(1 - 



7^ = (W) 2 



where and (+++) are those in Eqs. (107), (108) and (109). The matter energy is 

1 nr\ 1 d 



S uniform = _ Rg I 



In 

2 ve 



(!) 



■In 



7 S 



a 3 7s 



(203) 
(204) 

(205) 



Another special case is <p' = corresponding x' — » —00. In this limit, a = (3 and the q determinant Eq. (159) 
can be written in a simpler form 



detM; jn = 16a 2 cos 



jn 2 (™r\ + 4 sin 2 + 5 2 u;f + 8(5 2 + 25) sin 2 (%) 



4skr 



45 2 cos(^) 



(206) 
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and the q integral can be perform and the result of q energy contribution in the large M, N limit is 



£ v '=° = — V In 



2M 

KM/2 
,1/2 

<i / c?2; In 
Jo 



?7T 

16a cos — 

M 



KM/2 



4sin 4 (7rz) + 165 2 sin 4 (7rz/2) + 8(5 2 + 25) sin 2 (7rz) 
+ \J 16a 2 cos(ttz) + 4 sin 4 (7rz) + 165 2 sin 4 (7rz/2) + 8(5 2 + 25) sin 2 (7rz) 



where 

The matter energy in this case is 



4 sin 2 + 5 2 i^ 4 + 8(5 2 + 25) sin 2 (fa) 
165 2 cos(^) 



(207) 



(208) 



C' =0 = -|(£i|*'=o +6l^=o) + ^ V2 ^ In 



4sin 4 (7rz) + 165 2 sin 4 (7rx/2) + 8(5 2 + 25) sin 2 (7rz) 



+ ^ 165 2 cos(ttz) + 4sin 4 (7rz) + 165 2 sin 4 (7rz/2) + 8(5 2 + 25) sin 2 (7T2;) (209) 

D Compactification and dimensional Reduction 

A final noteworthy simplification occurs when one employs dimensional reduction on some of the q compo- 
nents by imposing exact Dirichlet conditions, q = 0, on the boundaries for those components. Then instead 
of the e trick, one can use the spin projectors to rearrange the bond structure in the way already used for 
the ghosts which do obey exact Dirichlet conditions. Let's rewrite the path integral in the case where this 
is done for all components of the g's: 

Tji = E /^^cxp|^(^+^- 1 -2^^- 1 )(lng+^ln^^)l 



exp 



exp 



exp 



exp 



exp 



-E^ 

1,3 

a x -< 



J 2 



i,3 " i,j 



R 



R 



(^ + i-^)(^ + i-^)--(^ +1 -qD 2 



(l-P[)(l-Pl +1 ) 

£(i - '■ >r - : l^ 1 p UiW + i 4)(4 + i - 4) + PiM'i 



i,3 



~ D 1 - p ^ -p?+i) + pUi- Pi-i)) 44 \ 

i,3 ) 



(210) 
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To understand this formula, focus on a solid line crossing n > 1 up spins. Then in the original formula there 
are n — 1 Gaussian approximated delta functions: 




(211) 



Integrating over n — 1 q's leaves a single q on the solid line to be integrated: 

f2^m\ d{n - 1)/2 f , / a f f a " J 

(— J /^(w?) E/^i---^ex P |-— g^j, ( 212 ) 

where we have discretized the single integral, <7 — > fe/i?, where the vector A; has integer components and i? 
is the corresponding compactification radius (note that one could have a different R for each component). 
The process of dimensional reduction is to keep only the contributions with all k = (R — > is one way to 
effect this), absorbing the factor {a/2-KmR 2 ) d ^ 2 in the coupling constant: 

*--(ssi?r* <213) 

More generally, one could stop before this last step and have a worldsheet formalism for the compactificd 
quantum field theory. Then one would need to explicitly implement the constraints that all k^'s on a solid 
line are equal and all k\ = on dotted lines (i.e. those that cross dotted lines). These constraints are easily 
imposed by introducing angular variables 6\ on each temporal link (j — Then the constraints follow 

by inserting 

n / 70j* ex p H • K(i - pi)+ i|A - - fcrl)] } • (2i4) 

ij 

Taking the limit R — > oo returns to the uncompactificd theory. The explicit infra-red cutoff R then serves the 
purpose of 1/e which need not be introduced. This gives an alternative infra-red cutoff in which conservation 
of transverse momentum is maintained throughout, and hence is superior (at least conceptually) to the e 
cutoff. Whether it is computationally superior remains to be seen. 
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